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Abstract 

In models of the Universe with extra dimensions gravity propagates in the whole space- 
time. Graviton production by matter on the brane is significant in the early hot Universe. In 
a model of 3-brane with matter embedded in 5D space-time conditions for gravitons emitted 
from the brane to the bulk to return back to the brane are found. For a given 5-momentum of 
graviton falling back to the brane the interval between the times of emission and return to the 
brane is calculated. A method to calculate contribution to the energy-momentum tensor from 
multiple graviton bouncings is developed. Explicit expressions for contributions to the energy- 
momentum tensor of gravitons which have made one, two and three bounces are obtained and 
their magnitudes are numerically calculated. These expressions are used to solve the evolution 
equation for dark radiation. A relation connecting reheating temperature and the scale of extra 
dimension is obtained. For the reheating temperature Tr ~ 10®GeU we estimate the scale of 
extra dimension /i to be of order 10“®GeU ~ 10“^cm). 


1 Introduction 

Brane-world scenarios with the observable Universe located on a 3-brane embedded in a higher¬ 
dimensional space-time have attracted considerable interest recently. Such models with matter on 
the brane can reproduce the main cosmological data [1, 2, 3, 4]. 

A general property of extra-dimensional models is that although ordinary matter is supposed 
to be conhned to a brane, gravity propagates in the whole space-time. This entails the effect that 
gravitons produced in reactions of particles on the brane can escape to the bulk. Graviton production 
is strong in the early hot Universe, and can alter the time evolution of matter on the brane and, in 
particular, the primordial nucleosynthesis. 

In this paper we calculate graviton production in a model of hve-dimensional Universe with one 
large extra dimension. Matter is supposed to be conhned to the 3D brane. Time evolution of matter 
in this model is described by the generalized Friedmann equation + 2pp + ■ ■ ■ [5, 6, 7]. We 

consider the period of early cosmology, in which the term quadratic in energy density is dominant, 
/i/p -C 1 (p is the normalized energy density on the brane dehned in (4), p = (—A/6)^'^^, and A is 
5D cosmological constant). 

Because the space-time is curved, a part of gravitons emitted in the bulk can return back to the 
brane [8, 10, 11] and bounce again to the bulk. In paper [12] an analytical method to show that a 
bounce is possible was developed. In the present paper we investigate further conditions of a bounce. 
Solving the combined system of equations of trajectories of the brane and of emitted graviton, we 
hnd conditions at which graviton can fall back to the brane. We derive an equation for the interval 
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of time between graviton emission and its return to the brane and develop a scheme to calculate 
times of returns of graviton to the brane for multiple bounces. We show that in the period of early 
cosmology the ratio of times to/^i of graviton emission to and its return to the brane ti to a good 
approximation can be expressed as a function of x = m(ti)/E(ti), where {E(ti), m(ti), p) are the 
components of graviton 5-momentum at the time ti when graviton returns to the brane. 

As an application of the above results, using the distribution function of emitted gravitons of 
paper [11], we calculate the the components of the energy-momentum tensor of bouncing gravitons 
is transverse to the brane, k is a number of a bounce). For the first three bounces we 
obtain explicit expressions for and estimate their numerical magnitudes. The expressions for 

the energy-momentum tensor of bouncing gravitons are used to solve the evolution equation of dark 
radiation [11, 13]. Solving this equation, we hud a relation connecting the reheating temperature of 
the Universe Tr and the scale of the extra dimension /x. Qualitative constraints on Tr and jj, are 
discussed. 

In Sect.2 we review two approaches to the 5D model. 

In Sect.3 we solve geodesic equations for gravitons propagating in the bulk. 

In Sect.4 we solve the combined system of equations for graviton and brane trajectories and hnd 
conditions for return of graviton to the brane. We calculate the interval of times between graviton 
emission and detection as a function of graviton momentum at the time of detection. 

In Sect. 5 consider multiple graviton bounces. We calculate the {nn) components of the energy- 
momentum tensor of gravitons falling to the brane. 

In Sect. 6 we present qualitative numerical analysis of the energy-momentum tensor and discuss 
solution of evolution equation for dark radiation. 


2 3-brane in 5D bulk 

We consider the 5D model with one 3D brane embedded in the bulk. Matter is conhned to the brane, 
gravity extends to the bulk. In the leading approximation we neglect graviton emission from the 
brane to the bulk. The action is taken in the form 



where X 4 = y is coordinate of the inhnite extra dimension, . 

The 5D model can be treated in two alternative approaches. In the hrst approach metric is 
non-static, and the brane is located at a hxed position in the extra dimension [5, 6]. We consider the 
class of metrics of the form 


dsl = g^^^^dx^dx^ = t)dt^ + a^{y, t)rjabdx^dx’^ + dy^ 


( 2 ) 


The brane is spatially flat and located at y = 0. The freedom of parametrization of t, allows to set 
n{0,t) = 0. The energy-momentum tensor of matter on the brane is taken in the form 


% = diag 5{y){-p,p,p,p}. (3) 

For the following it is convenient to introduce the normalized expressions for energy density, pressure 
and cosmological constant on the brane which all have the same dimensionality [GeU] 
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Reduction of the metric (2) to the brane is 


= dt^ + a^(0, t)r]abdx°'dx^. (5) 

The function a(t) = a{0,t) satishes the generalized Friedmann equation [6] 

H^{t) =+ {p+af + (6) 

where H(t) = d(t)/a(t) and Pw(t) is the Weyl radiation term [1], which below is set to zero. 

In the second approach the brane separates two static 5D AdS spaces attached to both sides of 
the brane. The metrics of the AdS spaces are solutions of the Einstein equations of the form 

J d2 

ds^ = -fi{R)df^ + —- + pjR^dx^dxa, (7) 

Ji\^) 

where 

f.(R) = 

Below we consider the case pi = p 2 and = 0. Trajectory of the moving brane in the R, T plane is 
given by parametric equations R = rb{t), T = Tb{t), where t is the proper time on the brane 

= - 1 , ( 8 ) 

dot is derivative over t. Reduction of the 5D metric to the brane is 


= —dt^ + rl{t)dx°'dxa- 

The function rb{t) satishes the generalized Friedmann equation [14, 15, 16] 


— + (p+ 


(9) 


( 10 ) 


Equations (10) and (6) with pwit) = 0 are of the same form, and a^(0,t) can be identihed with r^(t). 
Below we consider the case a = p [4], so that (10) takes a form 


=p^ + 2pp. 

rbj 

The normalized velocity vector of the brane and the normal vector to the brane are 

h 


= {v ^, n") = (ffe, h), = {n , n") = ± 


fin) 


, fin)fb 


Here 


Tb = e 


Vf + 

fin) 


where e = ±. 

In the following we choose with the sign 


n^ = - 


n 

fin 


, fia)n 


( 11 ) 


( 12 ) 


( 13 ) 



3 Geodesic equations in the picture with static metric 

Let A be parameter along a geodesic. Geodesic equations in the metric (7) are 


cPf 

(fx^ 

d?R 


rp dT dR _ 

^ ^^d\ d\ - ^ 


+ 2r 




= 0 


dx’’ dR 
^^~dX~dX 
'dRV . 


-^R 


df\ ^p, dx^-dx^ 
d\ ^ dA dA ” ’ 


( 14 ) 

( 15 ) 

( 16 ) 


We consider solutions of the geodesic equations even in A. Integrating the geodesic equations, one 
obtains [13] 


dT 




C° 


d\ pIX 


dR 

lx 


= iR^ {C^f + C'^'^ - 


-t2 


(17) 


where (G^, G“, C^) are integration parameters. 

Tangent vectors to a null geodesic satisfy the relation 


B 


9ab- 


dx"^ dx 
dX dX 


= 0 , 


(18) 


from which it follows that = 0. Eqs. (14)-(17) were solved with the initial condition that T(0) 
and i?(0) are located on the brane world sheet: T(0) = uito), R{0) = rb(to). Here to is the proper 
time of the point on the brane world sheet uito), Hito) at which the geodesic begins, i.e. the time 
of the graviton emission. 

The components of momentum of a graviton propagating along a null geodesic are proportional 
to the tangent vector to a null geodesic 


(p^, p^, p“) 










p2i?2 


(19) 


where cr = ±. Also we dehne er as G^ = e(r|G^|. Expanding the graviton momentum in the 
basis (n^(f), n^(t), ed(t)), where ed = 5^/pa, we have 


p^ = Ev^ + mn^ + p“ed. 


a 


where 


p^p^ = —+ p“ =0. 


The components of p"^ in the two bases are connected as 


j, EeJjRTlP -mH 

P = -^-5-’ 

pRb 


p" = rb{EEl — me^/pUlP), p“ = 


-R 


E = p^ e\/pmprb 


pRR 

pRb 


m = p^Hrb — p 


R 


P 

pn 

e^/jRTlP 


P^Tb 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


The components E and m depend on t through rb{t). 
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Introducing 




7 = 


and expressing 7 through E and m, we have 


CT-i 


= 1 - 




7 = 1 - 




.-R 2 


(/i2i?2)2pT2 


= 1 - 


^ Ei7-m6v//i2 + i72 y r^(t) 


R* 


(24) 


(25) 


If at a time t graviton is on the brane world sheet R = rh{t) {t is a time of emission, to, or a time of 
return of the graviton to the brane, ti, we take the basis {v^,n^,e^)(t) at the time t and obtain 7 
in a form 

^ ~ {Ee^/jETlP - mHY' 


(26) 


4 Bounce of massless particles in the period of early cos¬ 
mology 

We consider the radiation-dominated period of the early cosmology, when p/p 3> 1 or,equivalently, 
pt 1. Supposing that the energy loss from the brane to the bulk is sufficiently small to comply 
with the observational data, we neglect in the conservation equation for the energy-momentum tensor 
the energy flow in the bulk. In the period of early cosmology, in the model with extra dimension, 
from the expression for energy density of relativistic degrees of freedom 


P{T) = 


ft:^7rV(T)T^ 

180 


(27) 


{g*(T) is a total number of relativistic degrees of freedom [18, 19] ) it follows that 


p{t) g.{T)T^ 


rpA 


p{R) gAT,)Tt Tt 

For times ti and t in the region of early cosmology, from the Friedman equation one obtains 


(28) 


p{ti) V ^b{t) j t 

Following [11], with the use of ( 8 ), the equation for the brane trajectory can be written as 


(29) 


dn 


= e 


9 2 • 

P nn 


dn \fpRfTn 


2 2 


H 




(30) 


Integrating Eq. (30) with the boundary conditions rt = nito), n = n(to), we obtain the equation 
for trajectory of the brane 


p^{n{t)-n{to)) = e 


dr ^fJFTlP 


Rbpo) 


H 


= e 


nito) Tbit) 


+ 6 


dr 

'nito) ^ 


o/JrTIp 


H 


- 1 


(31) 



From the first integrals of the null geodesic equations (17) we obtain 


dl 


(32) 


Integrating Eq. (32) with the initial conditions R = ri,(to), T = Tb(to), we obtain the equation for 
for a null geodesic (graviton trajectory) 


- lY'dAT - Mk))- 


(33) 


If graviton returns to the brane at time ti, we have R = rb{ti). Combining Eqs. (31) and (33) and 
using Friedmann equation, + 2fj,p, we obtain an equation for rb(ti) 


[eeTe/j(l — 7 ) ~ 


( 1 . 

1 \ dr 

P + P ^ 

Knito) 

nih)) Jr,(to) 

VP^ + 2PP 


(34) 


Eq. (34) can be interpreted as an equation which determins the time of return of graviton to the 
brane ti for a given time of emission to. It is seen that Eq.(34) admits solution only if 


eeTe_R — (+) 

Expanding the integrand of (34) in powers of /r/p, we have 


(35) 


p + /i 


v/p2 + 2pp 


= 1 + 


1 /p 


2 Vp 




8 Vp 


+ ... 


where p(t) = p{to){j-b{to)/ri,{t)Y. Integrating Eq. (34), we obtain 


-^ 1 - 1/2 _ 


((1-7) 


nito 


1 ) 


p 


14 Vpo 


rb{to) rb{ti) 
rh{ti 


(36) 


nito 


-1 - 


_p 

11 Vpo 


rb{ti 

rh{R 
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-1 + 


The series in (36) is convergent. Introducing 


rbiti) 


to\ 

ij 


\ 1/4 


(37) 


and substituting po = piz where po = p(to) and pi = p(ti), we transform (36) to a form 


_ ^^- 1/2 _ = 1^(1-^^) 


[( 1 - 7 ) 


2\pJ 7(1-z) Vpi/ 11(1-^) 


+ 


(38) 


The function 


fk(z) 


z{l-z’^) 
k{l — z) 


is monotone increasing with the maximum at the point z = 1 equal to 1. 


(39) 



4.1 Conditions of the fall of graviton on the brane 

The sign of the graviton momentum component m{t) at the time ti at which graviton returns to the 
brane is opposite to that at the time of emission to. From (23) we have 

E r~.eT^/^PTJ^ (40) 
m ~ CT-ff — eR\/— 7 )^'^^. 

Condition (35) is satished in the following cases: 


(i) er — +, — +; e — +; (ii)er — e — +; 

(iii)er = -, e/? = +; e = -; (iv)eT = +, e/? = -; e = - 


Let us hnd in which case is realized one of the possibilities: either (a) mo = m(ro) >0, mi = 
m(ti) < 0 , or ( 6 ) mo < 0 , mi > 0 . 

• The case (i)(a). er = = e = + . 

E and m are 

E{t) ~ i/IPTJE - H{1 - 7)1/2 ^ ^ _ ^)i/ 2 . 

Conditions 

m(ro) Ho - ij- 7)^/2 > 0 , m{ti) Hi - iJhI + - 7)^/2 < 0 . 


are satished, if 


Ho 




gi 

TTVffF' 


Here Ho = H{to), Hi = H{ti). 

• The case (i)(b). = e/j = e = +. 

E and m are 


(41) 


E{t) ~ - H{1 - 7)1/2 ^ Q - 7)1/^. 


From conditions mo < 0, mi > 0 it follows that 

, < (1 - 7 )'/' < , 

Vh' + H^ VJhTh! 

or Ho < Hi, which is impossible, because to < E. 

• The case (ii)(a). ex = = —, e = + . 

Analogously to the case (i)(b) in the case (ii)(a) there are no solutions. 

• The case (ii)(b). ex = = —, e = + . 

E and m are 

E{t) ~ -^/JpTJE + H{1 - 7)1/2 ^ Q ~ -iL + /h^T7^(1 - 7)1/^. 

Solution mo <0, mi > 0 exists, provided (41) is valid. 
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• The case (iii)(a). = —, e/j = +. 

E and m are 

E{t) ~ + - H{1 - 7 ) 1/2 r^-H- + - 7)1/2 ^ Q_ 

Because m is negative and does not change its sign, there are no solutions. Analogously in the 
case (iv) m is always positive, and no solution exists. 

To conclude, we are left with the solutions of the types (i)(a) and (ii)(b), which are physically 
equivalent, because Eqs. (17) with er = = + transform to equations with er = cr = — under the 

change A —)■ —A. In the following we consider the case (i)(a). 

4.2 Relation between emission and detection times 

To solve Eq.(38) we need to transform the expression (1 — 7 )“i /2 —1 to a convenient form. Introducing 
X = E/m and using Friedman equation, Ef^ = + 2p,p, and (26), we have 


1/2 _ 1 = 

^/JpTJe-xH 

- 1 = 

1 + p/p- xa/I + 2p/p 


H -x^H-^ + p^ 


^Jl + 2p/p-x{l + p/p) 


- 1 


(42) 


Relation (42) is valid at the endpoints of graviton trajectory, at emission point and at points where 
graviton hits the brane. 

Because, as discussed in preceding subsection, xi < 0, the relation (42) written at time ti takes 
the form 

1 + p/Pi — a /1 + 2 /i/pi 


(1_7)-i/2_i = (i_|2.i|). ^_ 

\xi\{l + p / pi) + a/I + 2/i/pi 

In the period of early cosmology expression (43) approximately is 

pV2pI 


(43) 


(1-7) ^/^ - 1 ^ (1 - |a;i|) 


a/I + 2pl Pi + |a;i|(l + p/pi) 


(44) 


The advantage of this form of 7 is that we have extracted the factor {p/pi)'^. Now the Eq. (38) can 
be written as 


f hV 1-^11 

VPl/ (1 + kl|)(l + h/Pl) 




2/11(2:01) + ••• 


(45) 


or 


1 - kil 

(1 + kil) 


77 ( 2 : 01 ) 


—^ 2 / 11 (^ 01 ) + ... 
Pi/ 


(1 + p/pi), 


(46) 


Eqs. (38) dehne 2:01 — {to/tiY^* through the ratio Xi = mi/Ei and pi = p{ti), or, equivalently, 
the emission time to through the ’’mass” and ’’energy” of the graviton at the time of return of the 
graviton to the brane ti. In the following, for practical calculations, in the region p/pi < 1 we use a 
simplihed equation 


1 - |a^il 

1 + ki| 


- fiizoi). 


Domain of applicability and corrections to this equation are discussed in Sect. 6 . 
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(47) 



5 Multiple bounces 


To consider multiple reflections from the brane of bouncing gravitons we use matrix notations. 
Introducing 


and 


we have 


K = K-^ 


+ 1 -H/ii 

H/^i -./JPJJF+I 


P = P pR, 







( 48 ) 


( 49 ) 


( 50 ) 


The case with multiple bouncings is illustrated by the scheme 

/ Po \ ( pf = pT^ \ ( pf = pT^ \ 

Eo (to) ^ = Er' (E) ^ = e°^^ (ts) • (si) 

y rriQ j y rrPi = —j y j 

The left column corresponds to the emission time to- In fh® next brackets, in the left columns are 
momenta of ingoing particle, in the right ones are the outgoing. Under reflection from the brane 
momentum p“ parallel to the brane and energy E are conserved, transverse momentum to the brane 
m changes its sign: 

E ^ E, p^ ^ p'^, m —)■ —m. 

Momenta ~ C'^'^/ri,{t) are are rescaled when moving from one bracket to the next 


Pn^ _ n{tn-l) f tn-lV^^ 

pIlI “ U(tn) ^ Wn y* 


(52) 


Transformation from ”in” to ’’out” components within a bracket is given by 


Mout 


^out 

mOut 


1 0 
0 -1 


Etn 

m^in 


= LM^^. 


( 53 ) 


Let us consider the first two brackets in (51). Using (50) we express Eq and mo through Ei and mi 
and obtain 


Eq 

mo 


out 


= z^^KoKi 


El 

mi 


( 54 ) 


/ 


= 2 :, 


-1 

01 


fj2 

% + 1 

fjZ 


, 1 _ HqHi 

fj? ' 


Ho 




IJ2 

Esi + 1 


Ho 


4 + 1 


ih. 


4 + 1 


4 +1 


-Wf , 1 _ HqHi 

IE ~1E~ 


\ 

/ 



Q 



In the period of early cosmology, from the Friedman eqnation it follows that Ff//i ~ l/(4/rt). For 
small fit we can simplify the expressions in (54) as 








HoH, 1 f Ho HA 1 fh , to 
A ~ 2[h,^ Ho) ~ 2 \to^ h 


th 

fi 

Introducing 


we obtain 


/i 


Ht H, H^ ^ 1 [Ho HA ^ 1 

" ^ 2 + 2 [ Hi Ho) 2 


jj, V jj,^ 


h 

to 


A 


^01 + Ai 


-4 4 

^01 ^01 


= 


Z-^ ± 


Ho 

mo 


out 


= z, 


-1 f Ai 


01 


where 

For multiple bouncings we have 


Ai 

i^oi = i^^izoi) 


El 

mi 


MA = z^KoKiMr, 
= z^Ko KiLKi 


= Z^^lKoKiLKi 


where 


ZOn — Zoi Zi2 Z23 - ■ ■ Zn-l,m 

KoKiLKi...Kn-lKn = 


i’tn 

)”+Vi. 


V'o., 


Here 


i^On = ^^(WOn), 


±/ 


-1 


^On ^0,n—l Zn—lyn: 


v -1 


v -1 


v-1 w „-l. 


or explicitly 

%2 = Zq{ Zi2, Uo3 = Zoi Zi2 Z23, U04 = Zq{ Z12 Z2{Z34, ■ ■ ■ 

Fl'om (60) and (61) it follows that 'Uo, 2 fc+i < 1 Uoot > 1 ■ In the latter case 

^0,2fc(^0:2e) — ~IV'o,2fc(tto.2r)| — ~V'o,2t(tto,2ft)> 

and 


(55) 


(56) 


(57) 


(58) 

(59) 


(60) 

(61) 


(62) 


Eo = Zo,2kA+2kE2k - \ A,2kA2k) 
mo = Zo,2k{AAk\^‘^k - '4^0,2k^A- 



It should be noted that the functions Zn-i,n in processes with different number of bounces are different. 

(k) (k) 

If „ refers to the process with k — 1 bounces, different ^ are connected by the following 
relations 


,(fc) (fc-1) (fc-2) 

'k-l,k ~ %-2,fc-l ~ %-3,fc-2 — • • • 

— ^01 ; 

(63) 

,(fc) (fc-1) {k-2) 

'k-2,k-l ~ ^fc-3,fc-2 ~ %-4,fc-3 “ ‘ 

• • - 

~ ^01 ) 


(fc) _ (fc_i) 

12 ~ 




The distribution function of non-interacting gravitons in the bulk satisfies the Liouville equation 
without the collision term. If coordinates and momenta gravitons along a geodesic are parametrized 
by parameter A, i.e. f{x^{\),p"^{\)), we have 

/(fl(A„),p-’(A„)) = f(R(\j),p-*{M))- (64) 

In the case of the metric (7) relation (64) can be written as [11] 

f{fo,Ro,P^) = f{f,R,,p^R,/Ro), 

where p^ = (p^,p^,p). If the points (To,i?o) and {Ti,Ri) are on the brane world sheet, they are 
functions of the proper time on the brane. Temperature of the Universe, T{t), is defined through the 
proper time t via (28)-(29). 


5.1 First fall of gravitons to the brane 

We suppose that the distribution function of emitted gravitons depends on Eq = a/uIq -|- p^, mo 
and temperature Tq, i.e. = f°^^{Eo,mo,T(to)). The distribution function of gravitons emitted 
at time to and falling back for the first time on the brane at time ti is 

n^^\E,,m,,T) = r\Eo{E,,m,,Zoi),mo{E,,m,,Zo,),To{T,zoi)), (65) 


where 


Eo = Zoi^'PoiEi + 'iPoimi) ( 66 ) 

^0 = + i^oi^i). (67) 

In the period of early cosmology from (28) and (29) it follows that Ti/Tq ~ Pbito)/Pbiti) = Zoi- 
We obtain the distribution function at time ti as 

= f^\zQ^\Eip:+{zoi) + {zoi)), z^^^Eip:-{zoi) + mi'ip^{zoi)), z^^^T), ( 68 ) 

where 2:01 is determined as a function of Xi = mi/Ei by Eq.(45). For p/pi 1 the terms 0{p/pi) 

could be neglected and 2:01 is defined via (47). 

Condition that xo = mo/Eo > 0 takes the form 

^ ^~(^oi) - \xi\ip^{zoi) 

i)+{zoi) - \xi\p)-{zoi) 



Because the nominator of this ratio is positive, this condition is equivalent to tjj (zqi) — {zqi) > 

0, or 

l+_\^ 8 . 

To show that this inequality is satished, we wright 

l+\Xi\ X 1+klL l + l + klLs /fjnx 

5.2 Next falls of gravitons to the brane 

The distribution function of gravitons emitted at time to, which bounced off the brane at time ti 
and fall on the brane the second time at time t 2 is 

n^^\E2,m2,t2) = r''^^\Ei{E2,m2),-mi{E2,m2),U) = fm 2 ),mi(E2,m 2 ),ti) (70) 

= f°^\Eo{Ei{E2,m2),mi{E2,m2)), mo(-Ei(E2,m 2 ),mi(E2,m 2 )), to). 

Tracing the propagation of graviton, we obtain 


-1 ( 4 

V ^o"i 4 


-1 ( A ^-1 ( ^^2 ^ 


^01 


-^12 -^12 / V m 2 


— ( 2 : 012 : 12 ) 


-1 [ 'l /’02 t /’02 


%2 ^02 / V m 2 


where U 02 = 2 :oi 2 :;^ 2 ^ 

i’m = ^ [^^(2:oi)^^(^12) - ^■(^oi)^^(2:i2)] = ^ [^ 02 '^ ± ^ 02 ] • 
The time ti of the bounce is determined from Eq. (47) with z = Z 12 — (ti/t 2 )^^"^. 


1 - ^ 2 ! 

1 + \X2\ 


— frizu), 


where X 2 = (m2/i72)“. The time to is determined from the equation (47) with z = zqi — (to/ti)^'^"^ 


1 - kil 

1 + ki| 




where xi = (mi/Ei)®*^. Expressing xi through X 2 , we have 


Xi = x^ = 


ijj (2:12) +X2^+(2:12) _ ijj (2:12) - |x 2|^+(2 :i 2) 
^+( 2 : 12 ) +a;2^“(2:i2) ^+( 2 : 12 ) - k2|^"(2:i2)’ 


Substituting (74) in (73), we obtain 

f ( \= ^ ~ = ^+(^ 12 ) - k2|^~(^12) - ^"( 2 : 12 ) + k2|^^(2:i2) ^ 1 + ^ 2 ! 8 77C-X 

^ l + kil ^/’+(;2 i 2) - |a;2|^“(2:12)+^"(2:12) - |a;2|^+(2:12) l-|a;2|^^^‘ 


19 . 



Condition > 0 yields the constraint 'ijj {Z 12 ) — \x 2 \' 4 ^^{Z 12 ) > 0. This condition can be rewritten 
as (1 — 4)/(i + ^12) > 1^21, or equivalently, as 


(1 + \^ 2 \) ^ ^ 

“(l-k' 2|) ’ 

which is valid, because of (75). 

Condition mo > 0 is —(t/’o2 ~ k2|t/’o2) > 0) or 2:01(1 “ l^al) < Zi 2 i^ + \^ 2 \)) which is satished, 
because 

l + \x2\zf2 _ frizoi) 


1 + 

\ X2 \ 

^12 

1 - 

\ X2 \ 

000 


Using relations (72) and (75) we can show that 2:12 > ;2oi- 

The distribution function of gravitons (70) cab be expressed as 

f^‘^\E2,m2,T) = /°“*l°)(^02^(-E2l/>02 -"i2^(]2), ^cr2^ (^2^(]2 + "^ 2 ^ 02 ) , 2:o2^^)- 
In the case that gravitons have made two bounces, using (57) and (58), we have 


= Z0SK0K1LK1K2LK2KS 


-1 f 1^03 ' 4 ’o 3 


'^03 


Here 2:03 = 2:012:122:23, 1/’^ = and Mqs = 2:oi2:;^2^^23- The functions Zt-i^k are determined from 

the equations 

frizu) = ^ _ j2,gj^23> (’^S) 

Using the above relations, we can show that 2:23 > ;2i2 > 2:01. 

The distribution function f^^\m^,E^,T) is 

/(^^(ms, E3, T) = (2:o3^(E3i7o"3 + ^ 3 ^/;+^), ZQs\E3^/J+^ + rngi^os), ^o 3^T) (81) 

6 Numerical estimates and discussion 

We perform numerical estimates of the {nn) component of the energy-momentum tensor of incoming 
gravitons using the distribution function of paper [11], Qualitatively, the energy-momentum tensor 
of incoming gravitons at the registration time fi is formed by summing contributions from gravitons 
emitted at times to preceding the registration time ti. 

The distribution function of emitted gravitons is 


/(°^(m,p,to) = Bm^e 


210715’ 



where E = ^Jm? + p^, A is the weighted sum of relativistic degrees of freedom which contribute to 
the annihilation amplitude to gravitons [8]. The (nn) component of the energy-momentum tensor of 
emitted gravitons, TnZ^\to), is 

Tnnito) = j dmdp^/®(m,p,to)- (83) 

The (nn) component of the energy-momentum tensor of gravitons falling back to the brane is 

= j dmidpi^/(^^(mi,pi,ti), (84) 

where the distribution function of infalling gravitons is 

f”’(^)(mi,pi,fi,T) = 5mo(mi,Ei,^oi)exp ^njoi) 

Eo^rrii, El, zqi) and mo (mi, Ei, zqi) are dehned by (66)-(67). Substituting these expressions, we have 
= j ^ dmi ml ijEl - ml j dEiz^^ {Ei'ijj^i + mi^()~i)^exp | 

= 2 t^B j dmimj Ef - mj J dEiZQi{EiijQi - mi^/’w)^exp 

Introducing x = rrii/Ei, we express Tnn as 

= 27r5 J dxxVl -a;2 J dElEl{^l)Ql - x^p^lfzol%x) exp | ^^oi) | ^ 


First, we integrate over Ei in the limits (0, cxd), and below we consider integration taking into account 
lower and upper bounds. For Tin we have 


= 27,BT‘r(8) [' dxx^Vr^z^^ix A'H 

Jo yVoi—^Woi) 

For the energy-momentum tensor of gravitons which have made one bounce we obtain 

<iE2 j ^ + ™2V'o'2 


( 88 ) 


(89) 


(‘2i (“2^ (2'\ i2i 

where 2:02 = 2:12'^ Zqi ■ From the inequality zu > zqi it follows that t/’o2(2^12/^01 ) < 0- Instead, we use 
\' 4 ’o 2 \ = '4’02 {z^i /z^ 2 )- is expressed as 


For the energy-momentum tensor of gravitons which have made two bounces we have 

= 27^5 f dEs j dm 2 ml J El-ml ^03^(^3^703 + exp ^ ^ 

Jo J-E3 

rl 


(90) 


03 


T 


= 2 eBT^V{S) / dx x^Vl^z^i , 

Jo (m - ^^03)® 


( 91 ) 



where 2 ;o3 = 43 2^12 41- 

Integration over E in the integrals Tnn’^^^ is performed for E > T^in- 
Assuming that Tmin is in the region of early cosmology, i.e. 10 < p{Tmin)/l^-, taking 
and using (27), we hnd that 

min 2 ‘ 


~ 200 2 
(92) 


Using the relation ~ f^^pi which follows from the £t of cosmological data [4], we have ~ 
{pMpiY/S tt. For p = 10“^^ 10~^GeV condition (92) yields T^m ~ 10^ -i- lO^GeU. 

Because of the high power of E in the integrals for the main contribution to the integrals 

is produced from the region near the upper limit of integration Tmax- Provided T^m Fmax, we set 
Train = 0. The functions Zk-i,k{x) and the integrands for characteristic values of x are given in 
Table 1 and Fig. 1. For the following it is convenient to introduce the notations 



27r5T®F(8) / dxl^^\x) 


(93) 


Table 1: The functions 4-i k 


for different values of x. 


X3 = X 

.2 

.3 

.4 

.5 

.6 

.7 

.8 


.870 

.775 

.633 

.399 

.112 

.0034 

< 10-4 


.887 

.816 

.729 

.617 

.468 

.267 

.079 


.899 

.845 

.787 

.722 

.648 

.557 

.438 


^The ambiguity in and in A in (82) is due to incomplete knowledge of the contribution of dark matter. We 
assume that the mass of particles which form dark matter is in the interval (20 ^ 100)GeF. In the period of early 
cosmology these particles are relativistic. Phenomenologically the acceptable number of dark matter particles with the 
mass in the above interval is ~ 100 [19]. With the number of particle species in the non-supersymmetric Standard 
model g* ~ 100, the total number is ~ 200. Because of the high power of T estimates weakly depend on variations of 
this number. 





Figure 1: The integrands /c = 1, 2, 3 as the functions oi x = calculated for Tr/T 3> 1. 

The integrals were calculated numerically by substituting the values for Zij{x). For the 

integrals J dxl^^'^ we obtain 

[ dxJ(i) = 0.0415, [ dxl^^'^ = 0.0295, [ dxl^^'^ = 0.0023. (94) 

Jo Jo Jo 

To study dependence of z{x) on fi/pi, we calculate z{x) making use of (45) and taking into 
account the next order in p/pi- In Fig.2 the results for I^^\x) are compared for two values of 
p!P i = 0.1 and 0.01. Although the values of z{x) for p/pi = 0.1 and pjpi = 0.01 considerably 
differ, the values of integrals / dxl^^\p/pi = 0.1) = 0.0460 and / dxl^^\p/pi = 0.01) = 0.0415 are 
close. 

The limiting temperature at which the emission begins is the reheating temperature Tr. The 
emission energy Eq is bounded by Tr. Using the expression of Sect. 5 Eq = zonEipPtn ~ 
have 

Emax 'ErZOu/ iJ^On ^J^0n)i 

where ZQn{x) and = ip^{uon) are dehned in (59) and (60) correspondingly. For we obtain 

T-<■•)= 2 .BT^ (95) 

where 7 is incomplete gamma-function. 

If T/Tr < 1, the region producing the main contribution to the integral is 1 > Zon > T/Tr. The 
corresponding region ofa;is0<a;<O(l — T/Tr), where 1 — T/Tr 1. At small x the integrand 
decreases as a power of x, and contribution from this region is strongly suppressed. 

If T/Tr 1, the region producing the main contribution to the integral is 1 > zqu > 1 — T/Tr, 
where 1 - T/Tr < 1. 
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Figure 2: Dependence of on /i/pi, pi = p(ti), where ti is the time of the first return of 

graviton to the brane calculated using (46). The function calculated for p/pi = 0.01 (curve 

(1)) and p/pi = 0.1 (curve (2)). For p/pi < 0.01 the values of I^^\x) are very close to the case with 

p/pi = 0. 


The energy density of dark radiation satishes the evolution equation [9, 11, 13] 


^ + AHpo ^ --iTZ + TZ - TZ), 
at jj, 


where ^ 


TZ{T) = - 2 ^BT»r( 8 ) . 


T^™(T) = 27r5T®F(8) ■ 


8 7(8,T^/T)) 
105 F(8) 


In the period of early cosmology from the Friedmann and approximate conservation equations, 
p + AHp ~ 0, it follows that 

P_J__ 2/tVp,(T)T^ 

p Apt Ab{pMpiY 

Here we substituted M| ~ pM'^i [4]. Eq. (96) is transformed as 


dpoiT) A _ 2 

frpem I rpem rjiin\ 
dT rjnPD\^ ) ^nn)’ 


Explicitly we have 
dpD A 


dT T 


--pz. = 47rTT(8)H --+ 


TT 8 \ 7(8 ,Tzj/T) 


32 105 J F(8) 


— ^ j dxI^^\o. 


.7(8) ZokTu/T 


( 100 ) 
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Figure 3: The functions Zqi{x)Tr/T) calculated for Tr/T = 20,15,10,5 (curves 

from top to bottom) (plot (a)). The same for I^‘^'>{x)'-^{S^{zqiZi 2 ){x)Tii/T) (plot (b)) and 
/(3)(2,)^(8, i^zQiZi 2 Z 2 -i){x)Tii/T) (plot (c)). It is seen that as k increases the region of x, which produces 
the main contribution to is shifted to smaller x and larger z{x). 


Integrating (100) with the boundary condition pciTfi) = 0, we obtain 


8 \ 7(8,Tr/T0 
V 32 105 J r(8) 


For the ratio po/p 'Vfe have 


po = -T^ / dT'T 


,31X8)^ 

/i2%^ 




7 ( 8 , ZqiXrIT') 

r(8) 


( 101 ) 


Pd(T) 

m 


4725 A 


oTd r 


2^7r^g,{T){pMpi)^ J, 

4725 A 


7(8,Tr/T') 32 

-0.224 L —- + — > 

TT ^ 


''Tr/T 




y5 


F(8) 


-0.224 


dxl^^\ 


M8,ZokTn/r) 

F(8) 


7(8,1/) 

F(8) 


32 ^ 
H-/ 

TT ^ , 




T'^dT' 

( 102 ) 


Here we substituted M| ~ I^M^i [4]. 

For Tr/T > 20 the integral is practically constant and independent of T. The main contribution 
to the integral is produced by integration over the region of y near the lower limit. Taking Tr/T = 20 
and performing integration over p, we obtain 


32 

TT 




dy i{8,,y) 

F(8) 

7 ( 8 , zpky) 
F(8) 


~ 6.56 ■ 10"® 
~3.8-10“X 


(103) 

(104) 


For the ratio of energy density of dark radiation to energy density of matter we have 


Pd 

P 


3.5-10 


-5 


T; 


R 


{pMj 


pi) 




( 105 ) 




A typical order of constraint on magnitnde of the ratio p^/pin the period of early cosmology, which 
follows from primordial nncleosynthesis, is |pd/p| ^ 0.07 [20]. 

From the gravity experiments it follows that characteristic scale of extra dimension Vextr ~ 
is less than 10“^cm, or p > 2 ■ lQ~^‘^GeV [1]. For p ~ the estimate (105) gives Tr ~ 

2.3 ■ 10^ Gel/. This valne of Tr is signihcantly lower than usnally accepted Tr ~ 10^ -i- 10^ Gel/, 
indicating that large extra dimensions can appear with low reheating temperature. The estimate can 
be improved, if there is more complete cancellation between two terms in (102), or for larger values 
of p. Because of strong dependence of po/p on p, the possibility of larger p seems more plausible. 
For the reheating temperature Tr ~ lO^GeV the scale of extra dimension obtained from (105) is 
p ~ 2- 10“®Gel/ {p~^ ~ cm). For higher reheating temperatutes the value of p rapidly increases: 
p ~ T^. For larger p, in the integrals for Thn , the lower bound of integration Emin increases (see 
(92)) resulting in smaller magnitudes of the integrals T^. For Tr > lO^Gel/ this does not change 
the above results significantly, but for smaller Tr the effect of the lower bound of integration must 
be taken into account. 

The result of our calculations showing that account of gravity radiation leads to high mass scale 
of extra dimension can be attributed either to an insufficient accuracy of calculations (although our 
tests indicate stability of the result), or indicate that ’’too large” extra dimensions are incompatible 
with this class of models. 


7 Conclusion 

In this paper in a model with extra dimension we have calculated graviton emission to the extra 
dimension. Graviton emission is signihcant at the high-temperature period of the evolution of the 
Universe. The key point of the present paper is solution of the system of equations for the brane 
trajectory and for geodesic equation for graviton trajectory. For a given value of graviton 5-momentun 
at the time of graviton detection, we calculated the time of graviton emission. We obtained the 
recursion relations enabling, in principle, calculate the energy-momentum tensor of gravitons falling 
back to the brane which have made an arbitrary number of bounces. For the hrst three returns 
of graviton to the brane we obtained the explicit expressions for the energy-momentum tensor of 
the gravitons falling back to the brane and made their numerical estimates. Solving the evolution 
equation for the energy density of the dark radiation, we obtained a relation connecting the reheating 
temperature and the scale of extra dimension and estimated the scale of the extra dimension. 
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